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This  model  provides  a method  for  establishing  the  fragment  hazard  produced 
by  the  mass-detonation  of  stored  ammunition  stacks.  Fragmentation  characteris- 
tics used  as  input  are  derived  from  small-scale  arena  tests.  In  the  case  of 
projectiles,  the  small  scale  test  may  consist  of  one  or  more  pallets  positioned 
to  yield  a representative  sample  of  an  entire  stack.,. 

Hazardous  fragmentation  at  any  distance  is  currently  defined  as  having  a 
density  equal  to,  or  greater  than,  1/600  fragments  per  sqyare  foot  and  each 
fragment  making  up  the  density  having  a kinetic  energy  at  impact  equal  to,  or 
greater  than,  58  ft-lbs. 

The  unique  feature  of  the  model  lies  in  the  fact  that  a complex  trajectory 
Is  calculated  for  each  fragment  recovered  in  the  small  scale  arena  tests.  Using 
the  Dahlgren  main  computer,  approximately  200  fragments  per  minute  can  be 
processed. 

Past  tests  have  demonstrated  that  virtually  all  the  fragmentation  going 
down-range  Is  produced  by  the  ordnance  (projectiles,  bombs,  etc.)  on  the  face 
of  the  stack.  Fragmentation  from  the  ordnance  in  the  Interior  of  the  stack  is, 
for  the  most  part,  contained  within  the  stack.  When  a stack  is  detonated, 
fragment  jets  are  produced  between  adjacent  Items  on  the  face  of  the  stack. 

The  width  of  the  jet  is  dependent  on  the  method  of  stack  initiation.  When  all 
units  are  detonated  simultaneously,  the  jet  Is  typically  10  degrees  wide.  If 
only  one  or  two  donor  projectiles  are  initially  detonated,  the  jet  width  Is 
more  typically  20  degrees.  These  jets  are  referred  to  as  interaction  areas. 

The  greatest  densities  and  highest  velocities  are  produced  within  the  inter- 
action areas.  For  safety  reasons,  the  fragmentation  characteristics  of  the 
interaction  areas  are  used  for  Input  to  the  model.  The  Interaction  areas 
overlap  at  relatively  short  distances  down-range  and  can  therefore  be  added 
together  to  represent  the  cumulative  effect  of  large  stacks. 

Figure  1 shows  a single  pallet  of  projectiles  with  a one  degree  azimuthal 
slice  through  an  interaction  area.  The  one  degree  slice  has  been  selected  for 
mathematical  convenience.  The  slice  could  be  as  large  as  10  or  20  degrees 
depending  on  the  method  of  stack  Initiation.  From  arena  tests,  the  fragmenta- 
tion characteristics  can  be  established  for  the  full  90  degrees  of  ejection 
angle.  Individual  fragment  trajectories  can  then  be  calculated  to  establish 
hazardous  density  versus  range  In  terms  of  discrete  range  increments.  The 
hazardous  density  equation  is  shown  at  the  bottom  of  Figure  1. 
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Figure  2 shows  the  fragmentation  Input  to  the  program.  Each  fragnent 
recovered  In  the  arena  tests  Is  assigned  Its  own  specific  characteristics. 

The  ejection  zone  and  initial  velocity  for  each  fragment  Is  obtained  from  the 
instrumented  arena.  Fragment  weight  Is  measured  on  a scale  and  fragment 
average  presented  area  Is  measured  on  an  icosahedron  gage.  Subsonic  drag 
coefficient  (Mach  no.<.75)  has  been  correlated  with  the  ratio  of  maximum  to 
minimum  fragment  presented  area.  The  remainder  of  the  drag  curve  Is  approxi- 
mated from  the  historically  known  shape  of  fragment  drag  curves. 

The  program  has  two  options.  The  first  uses  average  values  for  the 
variables  shown  In  Figure  2.  The  second  Is  a Monte-Carlo  option  where  the 
uncertainty  in  ejection  angle  ( E) , initial  velocity  (V),  and  drag  coefficient 
( CQ)  are  simulated  by  sampling  from  appropriate  frequency  distributions. 

The  fragment  trajectory  for  each  fragment  is  calculated  using  a fourth 
order  Runge-Kutta  routine  (Figure  3).  The  trajectories  are  3-D  with  2-D 
wind  velocities.  Air  density  and  sound  speed  are  functions  of  altitude. 

The  sound  speed  Is  used  for  calculating  the  Cq  - Mach  no.  relationship.  The 
distance  and  kinetic  energy  at  impact  are  recorded  for  each  fragment.  The 
kinetic  energy  determines  whether  the  fragment  is  hazardous.  Summing  the 
number  of  hazardous  fragments  in  each  distance  increment  and  dividing  by 
the  area  of  the  one  degree  distance  Increment  yield  the  hazardous  density. 

This  is  compared  with  Hie  hazardous  density  criteria  to  see  if  the  criteria 
has  been  met. 

The  program  outputs  a table  as  shown  in  Figure  4.  In  addition  to  the 
number  of  hazardous  fragments  and  the  hazardous  fragment  density,  the  table 
includes  the  same  information  for  the  total  (hazardous  and  non-hazardous) 
fragments  and  the  ratio  of  hazardous  to  total  number  of  fragments. 

From  the  output  data,  hazardous  density  may  be  plotted  as  a function  of 
distance  as  shown  on  the  left  side  of  Figure  5.  In  this  example  the  number 
of  Interaction  areas  (NIA  - approximately  the  number  of  projectiles)  is  set 
to  50.  The  density  Is  proportional  to  the  number  of  Interaction  areas;  that 
is.  If  the  NIA  were  doubled  then  the  density  for  all  distances  would  also  be 
doubled.  Wind  has  a significant  effect.  A tailwind  will  Increase  both 
distance  and  Impact  kinetic  energy  as  shown  by  the  upward  and  right  shift  of 
the  peak  on  Figure  5.  The  approximate  Increase  In  distance  due  to  wind  is 
equal  to  the  wind  velocity  times  the  time  of  flight. 

The  density  versus  distance  data  may  be  translated  into  hazardous 
distance  versus  number  of  projectiles  as  shown  on  the  right  side  of  Figure  5. 
Knowing  the  hazardous  density  for  an  NIA  of  50,  the  NIA  for  a density  of 
1/600  fragments  per  square  foot  can  be  calculated  using  proportions.  The 
dotted  lines  on  the  plot  at  the  right  represent  a current  uncertainty  for 
very  short  distances.  Additional  work  with  small  ejection  angles  and  lighter 
fragnents  Is  necessary  to  evaluate  the  contribution  from  hitting  a standing 
man  at  short  distances. 
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The  essential  characteristics  of  the  program  are  as  follows: 

- Fragmentation  characteristics  derived  from  representative  small 
scale  arena  tests 

- Individual  3-D  fragment  trajectories 

- 2-D  wind  (horlzonal  plane) 

■*  4th  order  Runge-Kutta  Method 

- Average  value  and  Monte-Carlo  options 

- Air  density  a function  of  altitude 

- Sound  speed  (Mach  no.)  a function  of  altitude 

- Draa  coefficient  a function  of  the  maxiniimt  to  minimum  fragment 
presented  area  ratio 

- Handle  different  hazard  criteria 

- Output 

Hazardous  density  versus  distance 

Hazardous  distance  versus  number  of  projectiles  (bombs,  warheads, 
etc.)  on  the  face  of  the  stack  is  determined  from  the  density  tables. 


HAZARDOUS  FRAGMENT 

DENSITY 
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FIGURE  1 


RAW  FRAGMENT  DATA 


FIGURE  2 


FRAGMENT  TRAJECTORY 


FIGURE  3 


OUTPUT 


FIGURE  4 


HAZARDOUS 
QUANTITY  - DISTANCE 


